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Within a parton cascade approach we investigate the scaling of the differential elliptic flow V2(pr) 
with eccentricity e x and system size and its sensitivity to finite shear viscosity. We present calcula- 
tions for shear viscosity to entropy density ratio rj/s in the range from 1/4V up to 1/n, finding that 
the V2 saturation value varies by about a factor 2. Scaling of V2(pt)/cx is seen also for finite r//s 
which indicates that it does not prove a perfect hydrodynamical behavior, but is compatible with 
a plasma at finite r/fs. Introducing a suitable freeze-out condition, we see a significant reduction of 
V2(j>t) especially at intermediate pr and for more peripheral collisions. This causes a breaking of 
the scaling for both V2(pr) and the pr— averaged V2, while keeping the scaling of V2(pt)/ (v2)- This 
is in better agreement with the experimental observations and shows as a first indication that the 
rj/s should be significantly lower than the pQCD estimates. We finally point out the necessity to 
include the hadronization via coalescence for a definite evaluation of rj/s from intermediate pr data. 

PACS numbers: 25.75.-q, 25.75.Ld, 12.38Mh, 24.85. +p, 25.75.Nq 



The Relativistic Heavy Ion Collider (RHIC) has suc- 
cessfully shown that a transient state of matter at initial 
temperature T and energy density e well above the one 
expected at the phase transition (T c ~ 170 MeV and 
e c ~ 0.7 GeV/fm 3 ) has been created. In particular the 
large value of the elliptic flow v-i indicates that such a 
matter, called quark-gluon plasma (QGP), behaves like 
a nearly perfect fluid. In fact the dynamics of the bulk 
of plasma (i.e. for transverse momentum px < 1.5 GeV) 
is successfully described by ideal hydrodynamics [1], at 
least for the most central collisions [2|. At higher trans- 
verse momenta, due to incomplete equilibration, the hy- 
drodynamical behavior breaks down as confirmed by the 
saturation of the baryon to meson ratio and by the quark 
number scaling of elliptic flow i> 2 0, 0, S El ■ In this 
intermediate px region (1.5 < pt < 5 GeV) kinetic the- 
ory provides the most reliable approach and indeed par- 
ton cascade has successfully predicted the V2(j>t) satura- 
tion pattern for pt > 1.5 GeV [8|. Furthermore the cas- 
cade approach hints at a parton cross section significantly 
larger than estimated in perturbative QCD (pQCD) in 
general consistency with the observed nearly hydrody- 
namical behavior. On the other hand a minimum vis- 
cosity is imposed by quantum mechanical considerations 
[j| and more recently a study of supersymmetric gauge 
theory in infinite coupling limit [lfj has given a lower 
bound for the shear viscosity to entropy density ratio 
v/ s ^ 1/4-7T. All known substances, from water to a me- 
son gas, obey this bound and indeed all of them are sig- 
nificantly above it PJJ . A first recent evaluation of shear 
viscosity in lattice QCD (1QCD) is consistent with the 
lower bound H, 3] and show a mild evolution with tem- 
perature in the range of temperature covered by RHIC 
[ 1 31 ] . The description of the elliptic flow, the strong scat- 
tering of heavy quarks, the measurements of pt fluctu- 



ations, all point 11[ to a rj/s that should be close to 
the bound and much smaller than the one expected in a 
pQCD regime [3], (about 5-10 times the lower bound). 
A first attempt based on Knudsen number analysis of 
vij<i x has lead to estimate rj/s ~ 0.11 — 0.19 depending 
on initial conditions [l5j . Moreover it has been started an 
effort in developing viscous hydrodynamics 3, 13] which 
is indicating a significant reduction of v-i even for rj/s at 
the lower bound. The effect appears to be quite strong 
in the calculations of Ref. [I3| while is less pronounced 
in Ref. [l(| where data onpx averaged elliptic flow, (V2), 
are fairly well reproduced as a function of the collision 
centrality with rj/s ~ 0.1. On the other hand the low pt 
minimum bias V2(pt) seems to favor calculations with 
an even smaller rj/s. Similar findings for the dependence 
of the average elliptic flow {1)2) on shear viscosity have 
been reported also in the context of a parton cascade [r| . 
However a detailed investigation of viscous effects on dif- 
ferential elliptic flow V2{pt) within a transport theory is 
still pending. 

In this letter we present a study of the elliptic flow and 
its scaling properties as a function of pt at finite rj/s in 
the range (47r) _1 < rj/s < ir^ 1 . The analysis is based on 
a parton cascade approach and it is mainly focused on the 
intermediate px region, where kinetic theory automati- 
cally accounts for non equilibrium effects. The main idea 
is to keep the rj/s of the medium constant during the colli- 
sion dynamics. The parton cross section is rearranged ac- 
cording to the local density and momentum values. Sim- 
ulations have been carried out for a large range of impact 
parameters in Au + Au collisions at ^Jsnn = 200 GeV. 
Some simulations have been performed also for Cu + Cu 
for a first investigation of the system size dependence. A 
first issue that we discuss is the scaling of the V2 with 
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system size. We show that if rj/s is kept constant down 
to the thermal freeze-out (e ~ 0.2 GeV/fm 3 ) a parton 
cascade exhibits a v^/ e x scaling in the whole px range 
investigated (up to 3.5 GeV). Therefore the prediction of 
the scaling is not a unique feature of ideal hydrodynam- 
ics (lij . However experimentally the (in)dependence of 
v 2 / on the centrality of the collision and on the system 
size is indeed a delicate issue as raised by recent publica- 
tions from PHENIX [U and STAR [5lJ. We point out 
that once a suitable freeze-out condition is introduced at 
e c ~ 0.7 GeV/fm 3 a cascade approach at finite viscosity 
can account for the breaking of the scaling for v 2 (pt)/^x 
together with a persisting scaling for v 2 (pt)/ (^2), as ex- 
perimentally observed. 

A first attempt of our investigation has been also to 
put a reasonable constraint on the rj/s value of the 
RHIC fluid through the observed v 2 (pr) pattern for 1 
GeV< pt <3 GeV. A significant dependence of v 2 (pt) 
on shear viscosity is found with a reduction of the satura- 
tion value of nearly a factor 2 going from the lower bound 
to rj/s = 7T _1 . However a definitive evaluation of rj/s is 
entangled with the observation of quark number scaling 
in the same px range and hadronization by coalescence 
plus fragmentation has to be self-consistently included. 

The partonic transport approach at the present stage 
does not contain the different aspects of the dynamics 
and in particular it misses the effects of the fields, which 
have not yet been included in a partonic transport code. 
However it is certainly a powerful approach whenever a 
finite mean free path has to be considered and in par- 
ticular at intermediate px where the hydrodynamical be- 
havior breaks down. We have developed a 3 + 1 dimen- 
sional Montecarlo cascade for on-shell partons based on 
the stochastic interpretation of the transition rate. Such 
an interpretation is free from several unphysical draw- 
backs and particularly suitable for an extension to multi- 
particle collisions as pointed out by Z. Xu and C. Greiner 
[23J. The evolution of parton distribution function from 
initial conditions through elastic scatterings is followed 
by propagating particles along straight lines and sam- 
pling possible transitions in a certain volume and time 
interval according to the Boltzmann equation for two- 
body scatterings: 
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where J. = J. d 3 pj / (2ir) 3 2Ej, M. denotes the transition 
matrix for the elastic processes and fj are the particle 
distribution functions. 

For the numerical implementation, we discretize the 
space into cells small respect to the system size and we 
use such cells to calculate all the local quantities. In par- 
ticular we evaluate at each timestep the local collision 
probability and decide whether or not a collision can oc- 
cur by means of a Monte Carlo algorithm. We have per- 



formed several checks to test the validity of the code sim- 
ilarly to what thoroughly discussed in Ref. [23j , obtaining 
the same results. More specifically we have performed 
tests to choose a good discretization for convergency of 
the results for the elliptic flow that is the main observable 
analyzed in the present paper. The calculations shown 
are performed with cells of transverse area 0.5 fm 2 and a 
longitudinal size of Arj s ~ 0.1, where rj s is the space-time 
rapidity. Furthermore we have implemented the subdivi- 
sion (or test particles) technique which allows for a better 
mapping of the phase space. This is indeed necessary due 
to the smallness of the cell volume. Our tests have indi- 
cated that N — 6 test particles are sufficient to achieve 
stable results for the collision rate and the elliptic flow. 

In kinetic theory in ultra-relativistic conditions the 
shear viscosity can be expressed as 24 [ 



77 = 15 P ^ A 



(2) 



with p the parton density, A the mean free path and (p) 
the average momentum. Therefore considering that the 
entropy density for a massless gas is s = p{A — fi/T), ji 
being the fugacity, we get: 
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where <Jtr is the transport cross section, defined as 
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We use a pQCD inspired cross section with the infrared 
singularity regularized by Debye thermal mass tod 0] : 



da 
~dt 
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where s, t are the Mandelstam variables and mr> = 0.7 
GeV. 

Our approach is to artificially keep the viscosity of the 
medium constant during the dynamics of the collisions 
in a way similar to [2a ]. This is achieved by evaluating 
locally in space and time the strength of the cross section 
needed to keep the rj/s constant. From Eqs. ([2]) and (J3]) 
we see that assuming locally the thermal equilibrium this 
can be obtained evaluating in each cell the cross section 
according to: 
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with rj/s set from 1 to 4 in units of the minimum value. 

Partons are initially distributed according to a stan- 
dard mixture of the density of participant nuclcons (80%) 
and of binary collisions (20%) calculated with a standard 
Glauber model. The eccentricity of the system is there- 
fore similar to the one used in standard calculations. We 
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also start our simulation like in hydrodynamics at a time 
t = 0.6 fm assuming for partons with px < Po =2 GeV 
a thcrmalized spectrum and for p? > Po the spectrum 
of non-quenched minijets as calculated in j26[. In prin- 
ciple the transport approach allows for an investigation 
of the important issue of thermalization which should be 
strictly related to three-body scatterings [23| . Here look- 
ing at collective modes like the elliptic flow we implicitly 
assume that the results do not depend significantly on the 
details of the collision kinematics once the shear viscosity 
has been fixed. 

A first objective of our study is to investigate the scal- 
ing behavior of the elliptic flow with the initial eccentric- 
ity and the system size to see if such a scaling typical of 
a hydrodynamical behavior [3] persists also in a cascade 
approach. The interest for such a behavior is triggered 
by the recent observation by the PHENIX Collaboration 
[201 ] of a scaling of V2{pt)/ (^2) up to pr ~ 3 GeV, a 
region usually considered out of the range where hydro- 
dynamics should work. In order to allow for a compari- 
son with the results from ideal hydrodynamics of Ref. 
we have followed the evolution of the system considering 
a finite constant viscosity, as in hydrodynamical studies 
the zero mean free path condition is implicit for the en- 
tire evolution of the system. The hadronic re-scattering 
together with the formation and decay of the resonances 
are neglected. This is justified by the fact that the bulk of 
V2(pt) develops in the early stage of the reaction, i.e. well 
before hadronization sets in, as found by several theoret- 
ical approaches El [27L |28| and more recently confirmed 
experimentally [22j,l29(. In Fig. Q]the 1*2 As (left panel) 
and vi/k (V2) (right panel) in the central rapidity region 
(|y| < 0.35) are shown as a function of transverse mo- 
mentum pt for different impact parameters and the two 
systems Au+Au (filled symbols) and Cu+Cu (open sym- 
bols) at 200 AGeV with our cascade approach when shear 
viscosity is kept constant at 1/4tt. The dot-dashed and 
dashed lines are the results for Au+Au at b=7 fm with 
rj/s = l/2ir and rj/s = 1/ir respectively. The value of the 
constant in the right panel is set to k = 3.1, as in [20j ]. 
where k (V2) is assumed to be equivalent to the initial 
eccentricity e x . 

A first important result is the clear observation of the 
scaling both as a function of centrality and system size 
either for the V2{pt)/(-x and the V2(pt)/(v2}- In fact, 
since the parton cross section is re-normalized in order 
to keep a small constant shear viscosity, dynamical ef- 
fects related to the different density and temperature 
conditions that are reached at the different impact pa- 
rameters, are damped. This indicates that the scaling 
V2(pt)/ (V2), which is advocated as a signature of the 
hydrodynamical behavior [20| is a more general property 
that holds also at finite mean free path or shear viscosity 
at least for values close to the lower bound. Moreover the 
scaling is shown to persist also at higher px (~ 3 GeV) 
where not only the scaling but also the saturation shape 
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FIG. 1: 21 (left panel) and 
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(right panel) in the central 



rapidity region (|j/| < 0.35) for Au+Au (filled symbols) and 
Cu+Cu (open symbols) collisions at y/s = 200 AG eV . Dif- 
ferent symbols refer to cascade simulations at various impact 
parameters for rj/s = 1/47T. In the left panel also results for 
Au+Au at b=7 fm with rj/s = 1/2-7T (dot-dashed line) and 
rj/s = 1/n (dashed line) are shown. The constant in the right 
panel is set to k = 3.1 as in [20l |. 



is correctly reproduced by the parton cascade approach. 
Recently V2(pt) has been investigated also with viscous 
hydrodynamics [HI 17 1. It is interesting that a similar 
(but weaker) px dependence is found with a quantitative 
agreement with minimum bias data for rj/s ~ 0.1. This 
is in general agreement with our calculations, performed 
at various impact parameters, as we can see in Fig. Q] 
(right) comparing our results (symbols) with the shaded 
area (PHENIX data). We however are not aware of an 
explicit investigation of V2(pt)/£x and V2{pt)/ (1)2) scal- 
ing within hydrodynamics. Our simulations show a good 
sensitivity to the shear viscosity especially at intermedi- 
ate px where v 2 {pt)/^x drops of about 40% when increas- 
ing ?)/s by a factor 4 above the lower bound. While the 
^2/(^2) scaling was initially considered to stand for the 
V2/e x scaling and it is quoted as a further validation of 
ideal hydrodynamics [20(, latest results from STAR 21 1 
show that the V2{pt) scaled by the participant eccentric- 
ity e x is not independent on centrality. In fact the build 
up of a stronger collective motion in more central Au+Au 
collisions is observed. On the other hand, a good scaling 
with centrality is instead observed for V2I '(^2) ratio. This 
feature will be further clarified in the following. 



Effect of QGP freeze-out - We investigate the effect 
of a freeze-out condition on the elliptic flow. Freeze-out 
conditions are justified by the fact that, at a critical value 
for energy density, hadronization sets in and parton dy- 
namics is no longer acting. To take into account such 
an effect we stop the interactions among partons as the 
local energy density drops below 0.7 GeV/fm 3 , an in- 
termediate value in the range corresponding to a mixed 
quark- hadron phase [l[. Previous calculations were per- 
formed with a freeze-out condition at e = 0.2 GeV/fm 3 
which corresponds to the end of a mixed phase or roughly 
to an hadronic-thermal freeze-out. We have checked that 
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a freeze-out at 0.2 GeV/fm 3 is practically identical to 
consider no energy density freeze-out at all. This is in 
agreement with the observation that both theoretically 
and experimentally the elliptic flow does not develo p si g- 
nificantly during the hadronic stage P, [2^, 27, H, 2t| , 
hence we will refer to such a calculation as the one with- 
out fr( 
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FIG. 2: Differential elliptic flow for Au+Au at different im- 
pact parameters with (dashed curves) and without (solid 
curves) a freeze-out condition (e C rit = 0.7 GeV fm -3 ). 



When the freeze-out condition is implemented a size- 
able reduction for the elliptic flow is observed (see Fig. 
H]), especially for the most peripheral collisions and at 
larger p^- Correspondingly our results show that the 
scaling of elliptic flow with the initial spatial eccentricity 
is broken (see filled symbols in left panel of Fig. [3]) . In 
particular v%/e x varies of nearly 40 — 50% from b=3 fm 
to b=9 fm in the intermediate pr region (<3 GeV). The 
amount of such a spreading is consistent with the data 
reported by 2l| for the centrality selections — 10% and 
10 — 40%, with central collisions exhibiting a bigger el- 
liptic flow to eccentricity ratio than the peripheral ones. 
On the other hand, the scaling of v-ij (i^) with the im- 
pact parameter is still observed (see right panel in Fig. 
IS]). We are therefore driven to the conclusion that the 
breaking of the v^ipr)/^ scaling, as observed in Fig. [3J 
traces back to the freeze out physics, which deserves a 
deeper investigation. 

We are mainly focused on the pr dependence but it 
is of course interesting to look at the behavior of the 
averaged (vz). We however note that a comparison of 
the averaged vq. in our parton cascade with experimen- 
tal data should face two main limitations. One is the 
lack of finite quark mass effect that are known to reduce 
the value of Vi up to a pr of the order of the mass, an 
effect known in a hydrodynamical picture as mass or- 
dering. The other limitation is due to the absence of 
resonance formation and decay which are known to af- 
fect the elliptic flow especially for pions. Both effects are 
relevant at pr < 1 GeV and therefore should be taken 
into account to estimate the absolute value of (^2), as 
appropriately done in hydrodynamical models. Never- 
theless we are mainly interested in showing that a cas- 



cade approach is able to reproduce the observed trend 
of (1*2) l^x with the number of participant N part . Within 
the ideal hydrodynamics picture, the pr averaged elliptic 
flow is approximately proportional to the initial spatial 
eccentricity e x , leading to a centrality independent value 
of the (v2)/e x ratio. However recent measurements per- 
formed by STAR and PHOBOS [U H)j have pointed 
out significant deviations from the scaling. The trend 
experimentally observed (squares in FigHJ is recovered 
in our cascade approach at finite viscosity together with 
the scaling of 1*2 (pr) / (^2)- To compare with the experi- 
mental data we divide our results by a factor 1.2, which 
corresponds to fit the experimental value at N part ~ 170, 
see FigHl We mention that a breaking of the scaling for 
the average (1*2) is seen also without freeze-out condition, 
even if the effect of freeze-out reduces the absolute value 
and enhances the breaking. 
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FIG. 3: Same as FigfTJ filled symbols refer to calculations at 
r]/s — l/47r; open symbols are for b=7 fm (squares) and b=9 
fm (circles) calculations at 77/s = 1/n. The grey band in the 
left panel refers to the results of simulations without freeze 
out (see left panel in Fig. [T}. 
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FIG. 4: -^p- as a function of participant number for Au+Au 
collisions at y/s = 200 GeV in the central rapidity region 
\y\ < 1. Cascade results for 77/s = 1/47T are divided by 1.2. 
Squares are the corresponding data from [30( ] . 

As a last point we discuss how information on the 
viscosity of the RHIC plasma can be inferred from the 
V2{pt) absolute value. One has to consider that at inter- 
mediate pt there are several evidences for hadronization 
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via coalescence and it has been shown that due to a coa- 
lescence mechanism the parton V2 translates into a nearly 
doubled hadron v 2 [U [13] , therefore a definite evaluation 
of rj/s from i> 2 (jPr) data needs to include the coalescence 
plus fragmentation mechanism that should account for 
the baryon-meson quark number scaling. We refrain from 
using simple naive coalescence formula [33J here, con- 
sidering that it has been shown that space-momentum 
correlation and the freeze-out hypersurface can signifi- 
cantly affect the relation between quark and hadron V2 



34 , 1351 . 136j . It is therefore necessary a further develop- 
ment of the parton cascade approach that includes self- 
consistently the coalescence and fragmentation process. 
Nonetheless from Fig(5]we notice that for rj/s = ir^ 1 the 
parton elliptic flow with a quite small slope at low pj- sat- 
urates at about 6%. Even assuming a coalescence mech- 
anism in the hadronization phase, this value appears to 
be too low to reproduce the baryon and meson V2- This 
provides anyway an indication that a shear viscosity as 
high as 4 times the minimum value should be ruled out 
for the RHIC fluid and the viscosity is therefore quite 
smaller than pQCD estimates On the other hand 

the results with both 77/s = 1/Att and r]/s = l/2w could 
be quite close to the experimental data within a coales- 
cence picture. Such a range of values, to be narrowed in 
the next future, is sligthly larger than the first estimates 
with viscous hydrodynamics lg, 12 1 j slig htly below to the 
one based on Knudsen number analysis 15| and contains 
the best present evaluation in 1QCD 

PIT 



condition for the partonic dynamics is put at e ~ 0.7 
GeV/fm 3 then a sizeable breaking of the V2{pt)/^x scal- 
ing is seen while v^Ipt)/ (^2) still scales. This is in quali- 
tative agreement with the recent experimental data from 
STAR [2l| indicating that freeze-out of QGP dynamics 
with the consequent change of rj/s should be more thor- 
oughly investigated. As a final remark we notice that 
without any freeze-out condition the V2{pt) at parton 
level would be close to the data for rj/s = 1/4-7T, see 
Fig]?] On the other hand we consider such an agreement 
misleading because it would not be compatible with the 
enhancement of V2 due to coalescence and the observa- 
tion of quark number scaling 32]]. Instead the freeze- 
out condition seems to pave the way for a consistency 
among the different available observables on elliptic flow: 
the breaking of V2{pt)/^x, the persistence of V2(pt)/ (^2) 
scaling and the presence of a coalescence plus fragmen- 
tation hadronization mechanism acting at intermediate 
Pt- We therefore conclude that a safe evaluation of 
shear viscosity from the available data on V2(pt) neces- 
sitates a cascade approach that includes self-consistently 
hadronization by coalescence and fragmentation. Finally 
we mention that such an investigation can be strength- 
ened by a study of the fourth harmonic in the azimuthal 
anisotropy, i.e. the W4 = (cos(40)). A first analysis shows 
a stronger sensitivity to rj/s and especially a more critical 
dependence on the freeze-out dynamics respect to [37 1. 
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FIG. 5: Differential elliptic flow for Au+Au collisions at b=5 
f m i \y\ < 0.35 and rj/s = 1/4-7T (solid curve), rj/s — l/2rr 
(dot-dashed curve) and rj/s = 1/n (dashed curve). Results 
from cascade are compared with data from [20J (squares). 

Summary and Conclusions - We have investigated 
the dependence on the shear viscosity of the elliptic flow 
V2(pt) arid its scaling properties. As a first result we find 
that the approximate scaling of V2 {pt)I £x advocated as a 
signature of the perfect hydrodynamical behavior [20| can 
still hold also at finite viscosity and in a parton cascade 
approach. However such a scaling versus centrality and 
system size is present only if one makes the fireball evolve 
down to energy density e ~ 0.2 GeV/fm 3 corresponding 
typically to the end of a mixed phase. If a freeze-out 
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